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THE   EFFECT    CF    N0N-3PHERICi.L   SHA^E    ON   THE    MOTIOM    OF 

L   RISING  UlTDERVVATER   G/.S   BUBBLE 

(Preliminary  Memorandum) 

by  Max  Shiffman 

1  •      L^!^_rodu  ction 

The  problem  under  consideration  is  the  motion  of  a 
pulsating  gas  bubble,  produced  by  an  underwater  explosion,  as 
influenced  by  the  buoyant  force  due  to  gravitation.   This  pro- 
blem was  first  considered  by  C.  Herring,  using  a  perturbation 
method.   It  was  then  taken  up  by  G.  I.  Taylor,  who  integrated 
the  exact  equations  und^r  the  assumption  that  the  bubble  re- 
mains spherical  in  shape.   In  the  present  work,  this  assumption 
is  v/eakened  by  allowing  deviations  from  spherical  symmetry  but 
supposing  them  to  be  small.   But,  in  contrast  to  other  works  on 
non-spherical  shape,  we  do  not  assum.e  that  the  velocity  of  rise 
of  the  bubble  is  small. 

Experiment  shows  some  deviations  from  Trylor's  results. 
It  is  found  that  the  high  pressure  pulse  produced  v;hen  the 
bubble  is  contracting  and  near  its  minimum  size  is  larger  than 
that  predicted  by  Taylor's  theory.   Here,  qualitative  indica- 
tions are  given  that  non-sphericity  may  be  the  cause  of  these 
deviations.   Exact  statements  require  numerical  integration 

of  the  equations. 
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2.   The  mr-therncLticrl  fortnulr^tion  of  the  problem. 

The  surface  of  the  bubble,  which  is  rotntionr.lly  s-jn- 
metric  F.bout  v.   vertical  axis,  will  be  described  by  its  equation 
R  =  R(9)   in  polar  coordinates  with  a  specific  point  inside 
the  bubble  as  pole.   Determine  this  point  by  requiring  that  the 


'^xR   ) 


X^-^- 


r 


T\ 


expansion  of  R(9)   Legondre  polynomials  contain  no  ?,(cos  9) 
term,  i.e., 


.ZTT 


R(9)  cos  9  sin  9  d  9  =  0 


It  is  easy  to  show  that  such  a  point  exists,  and  is  in  fact 
unique  if  the  surface  does  not  deviate  markedly  from  the  spheri- 
cal shape.   The  equation  of  the  surface  is 


oo 


(1)     R  =  a(l  +  A  )     where 


X  =   <^  cx;  ?^(cos  9) 

•  •»):2.    J^  "  -^ 


Let  b  be  the  distance  of  the  pole  from  a  point  35  feet  above 
the  level  of  the  water.   The  quantities  a,  b,  cXg.o:^,  ,.  are  func- 
tions of   t   to  be  determined. 
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The  advmtr.ge  of  this  choice  of  the  polo  vill  r.ppear  :. 
later.   Here,  let  it  be  mentioned  that  if  one  retains  only  first 
powers  of  the  oc' s  and  neglects  product  terms,  the  volume  of  the 
bubble  is 
(2)  V=|Trr/^   , 

and  the  center  of  gravity-  of  the  bubble  is  rt  the  pole. 

If  a,  b,  oCr^,  . .  . ,  (X    , are  given  functions  of  t,  the 

hydrodynamical  state  of  the  water  at  any  instant   t   is  deter- 
mined by  a  velocity  potential  CD  (we  are  assuming  in  compressibil- 
ity of  the  v/ater)  which  satisfies  Laplace's  equation,  is  regular 
at  oo ,  r.nd  satisfies  the  boundary  condition; 

—  -3-- =  normal  component  of  velocity  of  the  surface  at 

the  surface  of  the  bubble.   Letting  x   be  the  angle  between  the 
radius  vector  and  the  normr.l  to  the  surface,  so  that 

the  bound;. ry  condition  becomes 

9cp         9  (p    . 

-  cosr  ^  •   +  sin  TT   —   '    =  R  cos  T  +  b  cosCC-r) 
or 
[^ 2_  +  tan  T   ■^--   =  R  +  b  cos  9  +  b  sin  9  tanT 

on  the  surface  of  the  bubble.   (The  dot  denotes  differentiation 

with  respect  to   t.) 
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In  tho  next  section  c.  nethod  of  successive  approxlmntions 
will  be  given  for  determining  the  potcnticl  function  sr.tisfy- 
ing  the  boundary  condition  (4).  Once  detornined,  the  notion  of 
the  gr.s  bubble  is  found  by  setting  the  pressure  of  the  water  at 
tho  surface  of  the  bubble  equal  to  the  gas  pressure.  The  pres- 
sure  ?   of  the  v/ater  is  given  by  the  Bernoulli i  equrtion 

where  ^  is  tho  density  of  the  water,  g   is  the  accelerr^tlon 
due  to  gravity,  and  the  extra  tern  b    '    occurs  because  of 

oz 

the  moving  coordinate  system  used  for  (D     . 

3.   The  detornin.'^ tion  of  the  velocity  potential. 

To  detomine  the  pot^ntirl  function  (P  satisfying  the 
boundary  condition  (4),  set 

CO 

(6)     <y  =  a<^,„.  iCf-^   *  iJ-'^n^'^n  ' 

where  'V    ,    Q -\  >    ^^o> >  ^    > ^"^^   potential  functions 

•     •  • 

not  depending  on  a,  b,  (a(X.  )  .   Substituting   (4),  rnd  using 
(3),  the  boundary  conditions  are 

^y  1+A  ^xe     '-^96 

-  -If^^T^  7F-1?-  =V-=«'     '"-^^ 

The  functions  (D    ,     CD,,....,   CD  ,....  can  be  found  by  success- 
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ive  approximj:.tion,  assuming  that  the  qurnti  ties  OCp,  OC,,  .  .  . ,  (V  ,  . . 
r.re  small  compared  to  1  and  retaining  successively  the  zero—, 
first  second,  etc,  order  term  in  the  oc's.   In  the  end,  only  the 
first  order  approximation  will  be  used. 

The  0—   order  approximation  for  W    ,     CP,,...   (X)  ,  .  .  . 
arc  easily  found,  and  we  can  sot 

n+2     P  (cos  0) 

where      l/J^>      ?/)-,>••••,  lif    »•••    ^^^   first   order   term.s.      Substi- 
tuting   (7),    and  retaining   only  first   order    terms,    one   obtains 


^  -r     la. 


(9)       )    ~  ~~^ /      =  3  A  cos   0     +5 sm   9 

-   -:i"i   =("+2)   P^(cos  0)  +  ^j  4-/^  sin  9  P^cos  9)n>2 
The  solutions  of  (9)  are 

Y-o    ^-  n+1  ^n       n+1 
/  ^   _  V_3n_  /  ^n-1   _  _njj._  )    g-    n 

The  expression  for  llj    ,  n  >  2,  is  too  complicated  to  set  down. 

But  note  in  (6)  that  (D      is  required  only  to  a  zero-  order 
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since  it  is  nultiDlied  by  {r.ix    )  ,  which  is  of  the  first  order 

n 

•  • 

compared  to  b.   This  follows  from  the  fact  thr  t  (aCX  )   is  sis;- 

n        ° 

nlfict'nt  only  when  the  bubble  is  contracting  rnd  near  its  mini- 

•  «     •  ■ 

mum  size.   During  this  str.ge  the  term  a  oc  is  (ci  OC   )    =   aoc +^0^11 
is  small  compared  to   b  because   a   is  small. 

Substituting   (8)  and  (6), 

00 

'sT   „n+2  ?  (cos  9) 


(11)  9=7   n^ n 

where 


o 


'1=^'  ■?"!  =1  -^"^2^ 

i.   =  -~    ( a  (X.  )  +  2  OC   a   +  -y~-   b  •^•^  -  y~^     n  >  2 
n   n+1       n       n       2      2n-l   2n+3    — 

The  potential  function  CD  can  be  obtained  to  a  second 

or  highter  order,  but  it  is  useless  to  do  so. 

4,   The  equations  of  motion. 

The  motion  of  the  bubble  is  nov;  obtrined  by  substitu- 
ting the  expression  for  CP  in  (5),  evaluating  at  the  surface  of 
the  bubble,  rnd  retaining  only  first  order  terms.   Exprnding  (5) 
into  a  series  of  Lcgendre  polynom.ials  and  equating  the  coeffici- 
ents to  zero,  the  result  is  a  series  of  equations  of  which  the 
first  is  simple  enough  to  write  down.   It  is 
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2    \' 

-  -  P-  b  +  4   b^   +  4-  a^  -  ^^   ' 

where  P  is  the  pressure  of  the  gas,  depending  on  a  according 
to  the  adlabatic  law.   This  cqur.tlon  does  not  depend  on  the  tx' s 
(to  the  first  order)  and  will  be  interpreted  in  the  next  section. 

The  remaining  equations  are  too  complicated  to  discuss. 
They become  simpler  if  we  suppose  that  the  equation  (l)  of  the 
surface  of  the  bubble  is  United  to 

(14)  X  =  ^2  ^2  ^^°^  ^^' 

i.e.,  ot.   =  o<   = =0.   In  that  case,  the  next  two  equa- 
tions arc  to  the  first  order, 

(15)  1^-^      +  2  g  a  (1  +  I  c(.^  )  -  I  a  b^^  =  0 

(16)  a^  ^2   ^   ^   -   ^'^*2  ^  ^2  ^-^  -  W     t)^)  ^2  "  "  1  ^^ 

The  three  equations  (13),  (15),  (16)  must  be  solved  numerically. 

It  is  convenient  to  set  down  here  the  expressions  for 
the  kinetic  energy  T   of  the  water.   Neglecting  quadratic  terms 
in  <u    ,  n  =  2,  3, ,  an  elementary  crlculation  yields 

(17)  T  =  2T^  Pi\^        a^  ^  i  ^^  ^^"i  *^2M  ^  quadratic  terms. 
(This  expression  for   T   can  not  be  used  to  derive  the  equations 
of  motion  in  the  Lagrange  form,  since  the  qurdratic  terms  are 
necessary  for  that.) 
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5.   Discussions  of  the  equations. 

In  this  section  we  shall  study  qualitatively  the  effects 
of  the  deviations  from  spherical  shape  on  the  upward  notion  of 
the  bubble . 
c)      The  sign  of  <3^^.   Experimental  evidence  and  some  theoretical 


graphs  of  Penney  shov/  that,  in  the  crucial  phase  when  the  bubble 

is  contracting  and  approaching  its  minimum  size,  the  quantities 
0^  2»  otp  become  negr. tive. 

b)  The  upward  velocity.   Writing  equrtion  (15)  as 

(^^b*)'  =  .  2  -  ^^  .  18      ^  ^    +   i   r  b  i 
a  "^ 

18  9    '  '  * 

the  two  terms  -=-  g  a  <>•  „  and   =  a  b  <x.  o  ^^e  positive  (b  is  nc- 

o  ti  o       d 

gative),  so  thet   (a   b)    is  not  so  large  negatively  (  and  may 

2 
a 

becom.e  positive  as  in  the  spherically  symmetric  one.   Thus,  the 

bubble  does  not  rise  as  rapidly  as  predicted  by  Taylor's  theory. 

1  *2 

c)  The  minimum  size.   If  one  groups  |  b   in  equation  (13) 

with  the  equation  (13)  can  bo  interpreted  as  follows:   in  addi- 
tion to  the  gas  pressure   P  opposing  the  contraction  of  the 
bubble  there  is  an  additional  opposition  represented  by  the  term 

7  v  b  .   Nov/  v;e  have  seen  that  b   is  smaller  than  that  predicted 

1   2 
in  the  spherical  case.   Therefore  the  opposition  |  <f  b   to  con- 
traction is  smaller  than  in  the  sphericr.l  case,  and  the  bubble 
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will  reach  a  smaller  minimum  size.   Thus,  the  effect  of  de- 
viations from  spherical  symmetry  is  to  cause  the  gas  bubble  to 
contract  to  a  smaller  size  than  that  predicted  by  Taylor's 
theory, 

d)   The  pressure  at  a  fixed  point  in  the  v;ater .   The  pressure 
at  a  fixed  point  in  the  v;ater  can  be  obtained  by  substituting 
(11),  (12),  into  (5).   The  result  is  an  expression  for   ?   in 
powers  of  -  .   We  are  Interested  in  the  pressure  when  the  bubble 
is  contracting  and  near  its  minimum  size,  i.e.,  when   a   is 
small,  while  the  point  at  which  the  pressure  is  derived  has  a 

relatively  large  distance  from  the  bubble.   The  doiiiinant  term 

a       "^ 
in  ?  is  therefore  the  first  power  of   -  ,  ©nd- *e-  have 

r  ' 


-3-   -  g  z   =  --^-i-   +  terms  -p-  , 


From  (13) 

,X8,  ^     .^     l^     .      l^.l^  .]      

v/here  P(a)  is  the  gas  pressure  inside  the  bubble,  and  hydro- 
static terms  have  been  ignored  in  (18). 

Nov;,  the  maximum  pressure  at  a  specified  point  in  the 
water  will  occur  at  or  very  near  to  the  time  when  the  bubble 

m 

has    its  minimum   size,    which  occurs  when     a      =     0; 

At  minimum  size 
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(19) 


V      )       o  4 


where   a   Is  the  ninimum  radius  of  the  bubble.   The  pressure  Is 

compounded  out  of  two  terms  -  the  internal  pressure  ?(a)  of  the 

1   *2 

gas  r.nd  the  dynrLinic  term  •:f  <?b  .   It  is  necessary  to  determine 

which  of  these  two  terms  is  the  most  important. 
The  energy  equation  is 


T   +   G(a)   +   I  TT  r.*^ 


(20)  T  +  G(a)   +   ^  ir  r.^  P   g  b  =  W 

where  Gr(r.)  is  the  internrl  energy  of  the  gr,s,  vV  is  c  constrnt, 

4  3 

end  ?  ^"   f  g  t)   is  the  v/ork  required  to  expand  the  bubble 

agrinst  the  hydrostatic  pressure   f*  g  b.   Neglecting  hydrostr.tic 

terms  r.nd  using  (17),  we  have 

(21)  t2  =  ¥-^-43-^ 


when   a  =  0. 


Substituting  (19),  and  using  G(a)  = 


{') 


4  r3 

3   ^- 


r- 1 


where  ?(a)  is  given  by  the  relation  ?(a)  =  ■--_■»- 


-or 


we   obtain 


(22)  ?   =   - 

r 


I  I'  F-(i  -fcc^    U^-  1)(1  -1^2) 


(") 


Using  equations  [10],  [12]  in  Taylor's  report  for  T.N.T., 

one  can  determine  the  value  of   a  which  will  maximize  the  co- 
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efflclcnt  of  -i   in  (22).   Tho  rosialt  for  the  ccso  <^  „  =  0  Is 


(23) 


a      =   1.30  M^/^ 


for  mnxin-um  possible  P,  where  M  is  the  mass  of  the  explosive. 
In  the  present  case  of  negative  o(„,  the  criticr.l,  value  of  a  is 
smeller  than  (23) . 

The  cr'. tical  value  of  a   is  to  be  compared  with  the 
minimum  possible  value  of   a  which  would  occur  if  there  were 
no  verticf-1  motion  of  the  bubble  and  which  is  determined  by 
G(a)  =  W.  This  value  of   a   is 

(24)  a  =1.25  H^/^      where  G(a)  =  W. 

The  minimum  radius  obta.incd  by  Taylor  in  his  calcula- 
tions is  3.3  tiie  (24),  or 

(25)  a  =  4.1  M^/^ 

A  comparison  of  (23),  (24),  (25)  shows  that  the  upward 
motion  of  the  bubble,  although  causing  the  contraction  to  take 
place  nearer  the  target  (if  the  target  is  situated  above  the 
explosion),  prevents  the  bubble  from  obtaining  a  much  smaller 
radius  "^nd  so  produces  r  correspondingly  weaker  pressure  pulse, 
The  most  imxportant  factor  in  obtaining  a  high  pressure  pulse  is 
to  rccach  as  small  c.   radius  as  possible. 

In  part  (c)  of  this  sect;' on,  we  showed  that  the  effect 
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of  devirtion   from   spherlcp.l    shr.pe   is    to  roach  r.    smr  ller   rtidius 
"c      than  otherwise.      Thus,    there  will  be   a  higher  pressure   at 
a    specified  point   in    the   water  then   that  predicted  by  Taylor's 
theory, 

6,      Conclusion, 

The   deviations   from  Taylor's   thoery  noted   In    section  5 
have   been  observed  experimentally.      Non-spherical    shape   of   the 
gas   bubble   gives   a   good  qualitative    explanation,    in    the   right 
direction,    of   these   observations.      To  ansvifer  whether   it  gives 
quantltrtive   agreement    with   the    obscrvrtions  would   require    the 
numerical    integration   of  equations    (13),    (15),    (16), 
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